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Abstract 



We carefully compute the gluon propagator in the background of a non-Abelian 



^ , Weizsacker-Williams field. This background field is generated by the valence quarks 

' in very large nuclei. We find contact terms in the small fluctuation equations of 

motion which induce corrections to a previously incorrect result for the gluon prop- 
agator in such a background field. The well known problem of the Hermiticity of 
certain operators in Light Cone gauge is resolved for the Weizsacker-Williams back- 
ground field. This is achieved by working in a gauge where singular terms in the 
equations of motion are absent and then gauge transforming the small fluctuation 
flelds to Light Cone gauge. 
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1 Introduction 



In a recent work, two of the authors have argued that parton distribution functions 
may be computed for very large nuclei at small x Q. They found that the valence 
quarks in these nuclei are the sources of non-Abelian Weizsacker-Williams fields. 
This work is similar in spirit to that of A. Mueller who considers the distribution 
functions generated by heavy quark charmonium-like bound states In refer- 
ence HI the general form of these fields was found. It was argued that to compute 
the distribution functions, one should compute the fields associated with these va- 
lence quarks and then color average over these fields treating the valence quarks as 
classical charges. The weight for this averaging was shown to be 

exp J d^xtp^ixt)^ (1) 

where p{xt) is the distribution of transverse valence quark color charge and /i^ is a 
parameter corresponding to the valence quark color charge per unit area. 
The current for the valence quarks was found to be 

Jli{x) = 5^'+6{x~)pa{xux+). (2) 

Here the Light Cone variables are 

y^ = {y^±y^)/V2. (3) 

The approximation which leads to this form of the current is that x « A~^/'^ so 
that the small x partons see a nucleus which is essentially Lorentz contracted to a 
thin sheet. This is the origin of the factor 5{x~) above. Further, since the small 
X partons are weakly coupled to the source, the source is "static" and only the + 
component of the current is important. 

The charge distribution in transverse coordinates is essentially random with a 
Gaussian weight as described above. This is true so long as we look at transverse 
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resolution scales which are large enough that there are many valence quarks per 
unit area, and small enough that it is much less than the confinement scale of 1 fm. 
In momentum space, this requires that we only consider transverse momenta which 
are in the range 

Kqcd « h « II ^ A^'^ . (4) 

The (Light Cone time) dependence of the charge density is a consequence of the 
extended current conservation law 

D^J^ = 0. (5) 

The current may be taken to be time independent if we fix the residual gauge 
freedom by the constraint that A~{x~ = 0) = 0. This point will be discussed 
further later. 

Under these conditions, it is easy to guess the form of the Weizsacker-Williams 
field. The charge distribution is confined to an infinitesimally thin sheet. Therefore, 
one expects that the fields on either side of the sheet might be a gauge transform 
of vacuum. We deduce therefore that Q 

= 
^- = 

A' = e{x~)a\xt) (6) 

where 

a^{xt) = -U{xt)VW\xt). (7) 
9 

The field a is a gauge transform of the vacuum. It generates zero magnetic 
field. There is a transverse electric field which is confined to the sheet x~ = 0. 
Gauss's law shows that this field solves 

V-a = gp. (8) 
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One can verify that the Light Cone Hamiltonian satisfies 

P- = (9) 

so that these classical fields have vanishing Light Cone energy. 

These fields are gauge transforms of vacuum fields with a discontinuity which is 
adjusted so as to reproduce the correct valence charge distribution. Because it is a 
pure gauge transform of the vacuum on either side of the infinitesimally thin sheet, 
we expect that the propagators are simple. In a previous paper, two of the authors 
argued that the fields which generate the propagator are just gauge transforms of 
the ordinary small fluctuation plane wave flelds Q. This argument is correct for 
matter fields such as scalars and fermions, but as we shall see is not quite correct 
for gluon fields. 

The problem with the gluon fields is that there can be contact terms in the 
gluon small fiuctuation equations of motion associated with the thin sheet of valence 
charge. These contact terms can generate contributions to the gluon field which are 
constant or grow linearly at large distances away from the sheet. The linearly 
growing terms will give rise to a constant electric field. In any case, neither of these 
terms were allowed for in the previous analysis which claimed to have determined 
the gluon propagator. 

The terms which were ignored in the previous analysis are very singular. A 
direct consequence of their bad behavior at infinity is the non-self-adjoint character 
of the Light Cone gauge equations of motion. 

The present work is organized as follows: in Section 2 we will see that the above 
mentioned terms are a consequence of the small fluctuation equations of motion 
in Light Cone gauge {A~^ = 0)- To resolve the ambiguities inherently associated 
with the singular behavior of the small fluctuation flelds, we choose to work in 
A" = gauge. In Section 3, we show that this gauge is non-singular and that 
the propagator can be computed directly. It turns out that the propagator in this 
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gauge is precisely that previously computed by McLerran and Venugopalan [Q. 
The propagator in Light Cone gauge is derived from the propagator in the non- 
singular gauge. Details of this derivation are the subject of Section 4. In section 5, 
we discuss the regularization prescriptions which enable us to generate a properly 
defined propagator in Light Cone gauge. We write down an explicit expression for 
this propagator. The propagator will be used in a sequel to this work to calculate 
the gluon distribution function to next order in a^. 
In Section 6, we summarize our results. 



2 The small fluctuation equations in Light Cone gauge 

The conventional approach to computing a Green's function is to write down the 
equations for the small fluctuation field in a fixed gauge. Denoting the small fluc- 
tuation field by 6(j), the generic small fluctuation equation is 

G-^6<Px = , (10) 

and the propagator is 

G(.,,)=E^^^^^^^. (11) 
To derive the equations for the small fluctuation field in Light Cone gauge {A~^ = 



0) we first look at equation (10) with A = 0. Expanding around our background 



field, we find that the equations for the gluon field become 

(^D{Afg^"' - L"^(A)D^(A)) 6A^ - 2Ff"'5A^ = 0. (12) 

In this equation, A is the background field which according to equation (|6|) is non- 
vanishing only for its transverse components. D^{A) is the covariant derivative with 
A as the gauge field (notice that = d^). We have assumed the gauge condition 
SA'^ (x^ = 0) = so that we can ignore any time evolution of the source of valence 
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quark charge. If this gauge condition were not chosen, then the extended current 
conservation condition would require that the source evolve in time, and this would 
generate a correction to the above equation. This background field generates no 
magnetic field, but it does generate a non-zero electric field, 



On either side of x~ = 0, the background field is a gauge transform of the 
vacuum. We might therefore think that the solutions to the above equation, in the 
two different regions of x~ , are simply different gauge transforms of the plane wave 
small fiuctuation field one would get in zero background field. Further, continuity 
at the origin could be assured by a matching condition. This is precisely what was 
done in the previous work of McLerran and Venugopalan |^] . 

Unfortunately, this procedure will not work here due to the presence of the 
background electric field term. Before seeing how this works, let us first express 
the small fluctuation equations in a slightly more transparent form (hereafter latin 
indices refer to transverse variables). 



D^{A)6A^' - D{A)^'D{A) ■ dA = -25^+5(x-)a^ x 6Ai - 26>'^6{x~)a^ x 6A-. 



The terms in the right-hand side of the above equation are called contact terms, they 
appear when the derivative d/dx~ acts on the classical field which is proportional 
to 9{x^). Notice that the last term vanishes in the gauge where 5A~ |^-=o= 0- 

The small fluctuation equation for the plus component of the above equations 
becomes 




(13) 



d+{d+dA- -D{Ay6Ai) 



26{x )q* X dAi. 



(14) 



This equation requires that 



d+6A- - D{Ay6Ai 



2e{x-)a' X 5A, 



,Q + /l(x+,Xt) 



(15) 
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where h is some as yet undetermined function. We see immediately that this imphes 
that the fields have singular behavior at infinity. We might remove the singular 
behavior at say x~ — oo but due to the step function, cannot remove it as 
x~ — > +00. In fact, the fields will behave like x~ as x~ goes to infinity. This 
singular behavior is sufficient to guarantee that the equations of motion are not 
self-adjoint on this space of functions. Further, the small fluctuation fields are 
discontinuous at the origin {x" = 0) and matching the fields there is more difficult. 
The singular behavior at infinity and the discontinuity of the fields at the origin 
were absent from the previous analysis of McLerran and Venugopalan Q . This was 
because they neglected to include the contact term at the origin which (as we have 
seen) is present in the equation of motion in Light Cone gauge. 

Due to this lack of self-adjointness on this space of functions, we have found no 
consistent way of directly integrating the equations of motion in Light Cone gauge. 
The most obvious problem is that the Green's function will not solve the constraint 
imposed by symmetry and by the Hermiticity of the underlying fields. Without a 
better way of defining the equations of motion, we are at a quandary and do not 
know how to properly construct the Green's function in Light Cone gauge. This 
problem with Light Cone gauge is well known in the relevant literature and various 
schemes have been adopted to circumvent it ||5|. 

However, and this is the core of the present analysis, we are compelled to 
compute the Green's function in Light cone gauge. There are several reasons for 
this: 

It is well known that only in Light Cone gauge does the parton model have the 
most direct and manifest physical realization Q]. Only in Light Cone gauge it is 
possible to construct a simple and intuitive Fock space basis on which our formalism 
relies heavily. Light Cone gauge (like any other axial gauge) is free of ghosts and 
negative-norm gauge boson states which makes the theory simpler [Q]. 

In addition, there are more advantages in choosing the A"*" = gauge. It can be 
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shown that in this gauge is the only large component of the fermionic current. 
The transverse and minus components of the current are proportional to and 
are therefore small [||. It is much easier to eliminate the dependent fermionic 
degrees of freedom in this gauge since we don't need to worry about inverting 
covariant derivatives. Finally, the gluon distribution function is most simply related 
to physical quantities measured in deep inelastic scattering experiments only in 
^+ = §. 

We will demonstrate in the following sections that the Green's functions in 
Light Cone gauge can be obtained by the following procedure: i) first, solve the small 
fluctuation equations in an alternative gauge {A~' = 0) and then ii) gauge transform 
the Green's functions thus obtained to Light Cone gauge. In the following section, 
we will show that the Green's function may be constructed in A~ = gauge and 
that the equations of motion for the Green's function are self-adjoint. In subsequent 
sections we discuss the mechanics of how the Green's function may then be gauge 
transformed to Light Cone gauge. When this transformation is carried through, 
there is in general a singularity of the transformation function. This singularity is 
regulated using the well known Leibbrandt-Mandelstam prescription We find 
that the Green's function in Light Cone gauge does not satisfy the constraints which 
arise from self-adjointness of the gluon field. This lack of self-adjointness is not 
irretrievably damaging and has been dealt with in the past for Green's functions of 
gluon fields in Light Cone gauge Q . Presumably the precise nature of the singularity 
is not important for gauge invariant quantities. When one chooses a prescription 
for regulating the singularity, if performed consistently in all quantities, it will give 
gauge invariant results for physical quantities. We are not aware that there is a 
proof of this, but believe it to be so. 
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3 The small fluctuations propagator in ^4 =0 gauge 

What is the virtue of writing down the equations of motion in A~ = gauge? In 
the small fluctuation equations for the transverse fields, the source term vanishes 
(the equation reduces to the sourceless Klein-Gordon field equation) 

2d+d- + D^{A)] SAi = 0. (16) 

This equation is obtained by appropriately fixing the residual gauge freedom in 
A" = gauge. One of the other equations is also sourceless, and is the Gauss law 
constraint which must be imposed on initial field configurations. This constraint 
commutes with the equations of motion. It is 

d- {d-SA+ - Di{A)6Ai) = 0. (17) 

We must now check that the remaining equation of motion is properly satisfied. 
If we use the Gauss law constraint that D^{A)SAn = 0, this equation is 

-2d+d- + D^{A)\ 6A+ = -26{x-)a' x SA\ (18) 

Now from the Gauss law constraint, 

SA+ = -^D{A) ■ SA. (19) 

Although SA^ is continuous at x~ = 0, the spatial covariant derivative operator is 
not. The discontinuity gives precisely the source term on the right hand side of the 
equation, and the equations are consistent. 

We can now solve the above equations of motion in the following way. We 
assume that the transverse field on the side where the background field vanishes 
is a simple plane wave. On the other side of the boundary where the background 
field is a gauge transform of the vacuum, we take the small fiuctuation field to be a 
linear superposition of plane waves. The linear superposition is done in a manner 
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to ensure that the fields are continuous across x~ = 0. This is possible because the 
equations above have no delta function singularity at = 0. 

The expression for the transverse components of the fluctuation field in A~ = 
gauge is 



A. 



0/3 



r/j0(-x~)T"^ + 0(x-) 



(2^)2 



X exp 



-I X 

2p- 



dhte-"^''' {u{xt)U\zt)TU{zt)U\xt)) 



al3 



(20) 



This solution is continuous at x~ = and was fully discussed in an earlier paper by 
McLerran and Venugopalan 

The above expression can be substituted in equation (pT[) to obtain the gluon 



propagator in A =0 gauge. There are four separate pieces for G'^^''^ ^ (x, y) which 



depend upon the signs of x and y . The overall expression for G is 
G?f^-'%,y) = -6., I ^:^[e{-x-)e{-y^)rfr^'^' 



(27r)4 p2 _ 



+ e{x-)e{y-)Ff{xt)F:'P\yt)+j J^d'zt 



(21) 



where F"^{xt) = yU {xt)TaU'^ (xt) j ■ It can be easily checked that this Green's 
function satisfies 

'-2d+d- + D^A)] Gkjix, y) = 5,, 5^^\x - y) 

as it must. This construction of the Green's function was previously performed by 
McLerran and Venugopalan in reference Q. However, they incorrectly interpreted 
it to be the result for A'^ = gauge, while we have shown here that it applies to 
A~ =Q gauge only. 
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For the pieces of the propagator which involve the component 5A^ we have 

< dA'dA+ >=< dA'-^D{A) ■ 6 A > (22) 

and 

< 6A+6A+ >=< ^D{A) ■ 6A -^D{A) -SAy- i6G++ (23) 

where in our convention < A^A^ >= —iG^v and < A^A^ > represents the corre- 
lation of fields in the path integral sense. 

4 On deriving the = propagator from the A~ = 
propagator 

In the previous section we derived an explicit expression for the small fluctuation 
Green's function in A~ = gauge. In this section, we discuss general arguments 
which allow us to compute expectation values of operators in one gauge in terms 
of the expectation values of operators in another gauge. The formal arguments 
will be illustrated by the specific example of the free gluon propagator. We are of 
course interested in these arguments because we want to find the gauge transform 
which will relate the Green's function in Light Cone gauge to the Green's function 
in A~ = gauge. 

Let 0{A) be some function of the field operators evaluated in Light Cone gauge, 
A'^ = 0. We will here consider QED and later generalize to the case at hand. We 
have to evaluate the expectation value 

J[dA] 5{A+) e'^0{A) 

We now introduce the identity in the form 

^ J[dA] 6iA- + d-A) 

J[dA] S{A- + d-A) ■ ^ ^ 
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Inserting this into the equation for 0{A), we find that 

J[dA][dA] 5iA-+d-A) 5{A+) e^^OjA) 
^ ' J[dA][dA] 6{A- + d-A) 6{A+) e^^ ' ^ ^ 

We now perform the gauge transformation 

A'' = A'' + d^A. (27) 

Since both the measure and the action are invariant under the gauge transformation, 
we see that 

< Q.^. I[dA][dA] 5iA-) 6iA+ - d+A) e^^OjA - dA) 

^ ' ![dA][dA]5{A-)6{A^ -d+A)eiS ' ^ ' 

Now doing the integral over A, we get 

([dA] 6{A-) e^^O{A-di^)) 

< 0{A) >= ^ ■ (29) 

^ ^ J[dA]5iA-)e'S ^ ^ 

Before proceeding any further, we will be more explicit and illustrate the above 
approach by an example: we show that the free propagator in the absence of a 
background field in = gauge can be obtained as a gauge transformation of the 
same quantity in A~ = gauge. 

Let us first review the construction of the propagator in A" = gauge. Hence- 
forth all fields in this gauge will be denoted by a ^ to distinguish them from the 
corresponding fields in Light Cone gauge. The action for the vector field is 

S = ^{d-SA+ - V ■ SAf + ^SA'SijOSA^ , (30) 

where the operator □ is the D'Alambertian in Light Cone coordinates. Prom the 
above expression for the action, we distinguish two classes of small fluctuation fields 
that satisfy the equations of motion. We can see this from the following argument. 
To find the solution that minimizes the action let us perform the shift 

6A+ ^ SA+ + • 6 A , (31) 
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to obtain a purely quadratic action in 6A~^ and 6Ai. The minimum thus corresponds 
to the fields that make the action vanish. We label the fields as follows: 
type I fields satisfy 

d-6A^-V-SA = 

aSA' = (32) 

whereas type II fields satisfy 

SA' = 

d-SA+ = 0. (33) 

The contributions of different components to the free propagator are then 

. . . . S'^ 

< SA'SA^ > = i— 

< dA'SA+ > = i 



The second term in the parenthesis above corresponds to the contribution of the 
expectation value < SAfj6Afj > to < 5A~^5A'^ >. Notice that the contribution 
from terms like < SAfSAfj > vanish. If we now combine all the above terms 
together we obtain the familiar expression for the gluon propagator 

where n*'^ = S'^'^. 

Now perform the gauge transformation 

5 A" = SAi" + d^A . (36) 

Since we are gauge transforming to Light Cone gauge, we require 

d+A = -SA+ . (37) 
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We therefore obtain 



5A+ = 

5 A- = -^5A+ 

6A' = 6A'-—6A+. (38) 

We now use equation ( p9| ) and the previous expression for the propagator in A~ = 
gauge to derive 

G>». ^""^') (39) 

rC \ rZ ' Ti / 

where = 5^^ . This result is the well known expression for the free gluon propa- 
gator in Light Cone gauge. 

For the case at hand, small fluctuations around a background field, the analysis 
goes through almost exactly as before. There are two essential differences. The 
first is that for the type of background field we consider, the gauge transformation 
property of the transverse field is 

5A' = 6A' + D'{A)A. (40) 

One should note that the particular form of the gauge transform above is valid 
only because we are interested in small fluctuations. We are throwing away terms 
which correspond to higher powers of the small fluctuation field. For the plus and 
minus components, the transformation involves ordinary derivatives. We see that 
the analysis above goes through in the same way it did for QED. 

The other difference from the QED case is that the structure of the small 
fluctuations analysis is more complicated. There are modes (labeled type I) which 
involve both 6A^ and 5A'^ as before but these modes now couple to the background 
field (see equation (|3^ ) with 5* replaced by D''{A)). There are also modes (type II) 
which involve only 6A~^ and therefore do not couple to the background field since 
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SA^ = 0. The contribution of the type I modes to the propagator arises from the 
distortion of the smah fluctuation fields by the background field. The contribution 
of the type II modes to the gluon propagator in momentum space is 

1 



6G 



(^^) " {k- 



\2 



- ^(^-2'-] (41) 

The second of the two expressions in the above equations is algebraically equivalent 
to the first but corresponds to a different way of regulating the singular behavior of 
the type II modes at = 0. In the next section, we will see that the latter choice 
is the one aptly suited to regularization. 



5 Regularization prescriptions and results 

We shall now use the method described in the previous section to derive the gluon 
propagator in Light Cone gauge from the propagator in A~ = gauge (equa- 
tion (^l|)). Note that the gluon propagator here is the gluon propagator in the 
background field. It will shortly become clear that we need to choose an appropri- 
ate regularization prescription to obtain the final expression for the propagator. 

If we have a fluctuation field in = gauge , the fields in j4+ = gauge are 
given in terms of the tilde fields as 

6A+ = (42) 
SA- = -^SA+ (43) 

6A' = 6A' - D\A)-^6A+ . (44) 

To derive these equations, we note that a small gauge transformation with a fixed 
background field A, where A has only transverse components, is 

6A+ + (45) 
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5 A' 5 A- + d~K (46) 

5A' 5A' + D{AYk. (47) 

We have used the fact that in A" = gauge 

d-5A+ = D{A)i5Ai. (48) 



The transverse fluctuation field in Light Cone gauge is then equation (^4|) which 
we rewrite in matrix component notation as 

5Af{x) = 



0-04- 



5Af{x). (49) 



To see how the covariant derivative acts on the field, notice that Dk5Ai:{x) = 
dkSAk — ig[Ak,SAk] where Ak{x) = 6{x~)^U{xt)dkU^ (xt) is the background field. 
Then the commutator of the background field with 5Ak will exactly cancel the 
derivative acting on U's. The result is that the covariant derivative will reduce to a 
simple derivative which acts only on the exponential part of the field. 

Following the discussion leading to equation (^), we obtain an implicit expres- 
sion for the Green's function in Light Cone gauge which we can write as the sum of 
the contributions of the type I and type II modes discussed in the previous section. 
The former can be written as 

(x) \ ^-^+ y (y) 



/IN ^7 / 1 \ 'I 



where Glff^j-^ ^ {x, y) is the Green's function in A" = gauge given by equation (21). 
The subscript / in the above equation denotes the contribution of the type I class 
of small fluctuation flelds. 

We must add to this expression the contribution from the gauge transform of 
(/;■) ■ This can be expressed as 

16 



where 

(*<5-);«V.',,.,) = (.J^*(.JV^7 l^^l^ . (52) 

It is now straightforward to construct similar expressions for non-diagonal com- 
ponents of the Green's function. Even though only transverse components of the 
Light Cone gauge Green's function are needed to compute the gluon distribution 
function, we will give an implicit expression for all components for the sake of com- 
pleteness. It is 



d. 


-d+ 




1 




-d+ 




1 




-d+ 


A 





(X) \ d-d+ J (y) 



{y) 

&^{x,y). (53) 



the Green's function in A~ = gauge. From here on, we will concentrate on the 
transverse components only. 

It is evident from the above equations that in relating the fields in the two 
gauges, one has to define and These operators are not well defined 

unless some prescription is adopted for treating the singularity at /c"*" = and at 
k" = 0. Presumably, there is no physics in these singularities -when one computes 
gauge invariant quantities, different prescriptions should give the same results. It 
has however proven most convenient to choose a Feynman like ie prescription for 
regularizing this singularity, the so called Leibbrandt-Mandelstam prescription [^. 
This prescription however spoils the Hermiticity of the field 5A^ . This is not a new 
problem since even in the absence of our background field, such a prescription will 
destroy the reality properties of A~ in = gauge. It is indeed unpleasant, but 
it does make for a Green's function which is manifestly causal. The price to be paid 
for a small lack of Hermiticity of the fields is small compared to the price of dealing 
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with non-causal Green's functions, albeit all such ambiguities must come to naught 
when computing gauge invariant quantities. 

We must now establish our prescriptions for defining the inversion of d^. We 
first begin with . For this operator we use the ordinary Leibbrandt-Mandelstam 
prescription 

k- ^ k- +ie/k+' ^ 
For the operator 1/5"^, we choose a prescription which guarantees that SA" \^-^q= 
0. This integration choice is in fact a further specification of the gauge and fixes the 
residual gauge freedom completely. This is because if we work in A'^ = gauge, we 
can still make x~ independent gauge transformations. We can therefore adjust the 
value of A~ at any fixed value of x~ which we choose to be x~ = 0. The virtue of 
this particular gauge choice is that it makes the source of charge time independent, 
since in this case the extended current conservation law is 

d-J+ = 0. (55) 

Our prescription will therefore be 

■^F{x-) = £ dx'-F{x'-). (56) 
Let us see what this corresponds to in momentum space. If we define 

F{k+) = J dx-f^^^""' F{x-) (57) 



then 



r 

JO 



dx'-F(x'-) = / — — '— (e-^'^^^ -1)F(A;+). (58) 



We see that with this integration prescription, there is never any singularity in 
Fourier space as A;"^ 0. 

The reason why this prescription is not often employed is because the types of 
fields one ordinarily considers are oscillatory at infinity, and therefore vanish in the 
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Riemann-Lebesgue sense. In such cases it is more convenient to choose = 
or = — since in the sense described above one can ignore the contribution 
from infinity |jlO|. In the Leibbrandt-Mandelstam prescription, positive frequencies 
propagate forward in time and negative frequencies propagate backward in time. 
For our problem where we have a nucleus acting as a source of charge, the fields do 
not vanish at x~ = ±00. There is therefore no advantage to choosing an integration 
prescription with one of its integration limits at either ±00. 

We can now combine terms together to get an expression for the propagator. 



The derivation of Gj,(//) from equation (51) is rather subtle and we shall elaborate 



a little on it. We can easily do the and kt integrations in equation ( |5^ to get 

dk~ e-*^~(^"^~J'^) 



6Gtr\Jx,y) = TaTa6'-^\xt - yt)S{x - y ) 



It is clear that we have locality in both transverse and minus coordinates. Hence 
there will not be any contributions from the regions 9{x~)9{—y~) and 9{—x^)9{y^). 
Then 



(D^)(..)(Z?^-)(,)<5G++ 



(-X )9{-y )TaTa didj5^'^\xt- yt) 
+ 9{x-)9{y-) Fa{xt)Fa{yt) d,d,5^^\xt - yt) 

_ f dk^ ^-'ik-{x+-y+) 

We now write the ^-functions in their Fourier representations and replace 



1 



1 / 



k+_\ 

{k-y ' k^ \ {k-y '^IFj ■ 



(59) 



(60) 



Finally, making use of the identity 



dk^ 



+ \ -ik+{x -y ) 



■{2k+)e 



d f dk 



dx 



27r 



ik'^{x —y ) 



d r dk 



dy- 



2-K 



^-ik+(x -y ) 



(61) 
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we can write equation (^) as 



G 



(27r)4 p2 



le 



^2 



PiPj 



P P 



l) + ^!^(e-ip+?/-_l) 



I PtPiPj ^^ip+x- _ -j^^j^g- 



1] 



(62) 



The derivation of an analogous expression for Giji^j-^ in equation (|5^) follows 
directly from the prescription rules discussed above. The reader may confirm that 
each term in equation ( p^ ) cancels an identical term in Gjj(/). Our final expression 
for the Green's function in Light Cone gauge is then 
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where = p'^ -\ — ^p-^ ■ Equation (p^) is the main result of this paper. In a 
forthcoming work, it will be used to compute the gluon distribution function in the 
Weizsacker-Williams background field to O(a^). There we will show explicitly that 
when computing the gluon distribution function from the expression (^), the only 
prescription needed to deal with poles is the Feynman (causual) ie prescription. 
In other words, the nature of the additional p^ poles (besides the usual one in 
p'^ — ie) in the Green's function is unimportant. Keep in mind that there are no 
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additional poles in since their residues vanish. As a result, we never have to use 
any additional ie prescription to compute the gluon distribution function to this 
order. We are not sure if this happens in all orders of perturbation theory or if it is 
an amusing coincidence specific to 0{a1). 

6 Summary 

In this paper we have computed the propagator for transverse gluons in the Light 
Cone gauge = and in the presence of a background Weizsacker-Williams 
field produced by the valence quarks of a heavy nucleus. We have argued that 
such a propagator can be obtained from the propagator in the A~ = gauge by a 
gauge transformation. This method has the advantage of avoiding the difficulties 
introduced by the singular behavior of the fields at x~ = ±oo caused by the presence 
of the sources. The gauge transformation requires a choice for the inverse of the 
operators d~, 5+. We chose these to be such that the gluon fields satisfy the 
boundary condition SA~\^-^q = and the usual Feynman causal behavior. The 
expression for the gluon propagator in = gauge is necessary if we want to 
extract from it the information about the gluon density which can be compared to 
experimental measurements. 

In the above analysis, there are still many subtle questions remaining about 
the independence of physical quantities from the precise choice of the regularization 
prescription for the small k~ and A;"*" singularities. We have shown that there is 
a gauge where there are no singularities and that the subsequent singularities in 
the Light Cone gauge Green's function arise from the gauge transformation rela- 
tions. Then, it is only sensible that the physical quantities computed here will be 
independent of the nature of these singularities. 

At the very least it seems we have a consistent prescription. It should be noted 
that because of the gauge choice, there is no ambiguity about the limit /s+ 0. The 
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propagator is smooth in this hmit as a consequence of the way we have specified the 
residual gauge degree of freedom. If it were only the free transverse propagator which 
was of interest, we would have a non-singular propagator at small k'^ . The problem 
we have with the distribution function is however at small . These singularities 
have been regularized in the only causal way possible, by an ie prescription. The 
question of physical interest is whether these small A;~, or correspondingly, large fe"*" 
singularities cause any problems at higher orders. It has been argued extensively 
in the literature that the prescription we chose for k~ is essential to prove the 



renormalizability of Yang-Mills theories in Light Cone gauge 1 11 ] . 
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